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We present a new method of determining the anisotropy of the gap function 
in layered high-Tc superconductors. Careful inelastic neutron scattering measure- 
ments at low temperature of the phonon dispersion curves in the (100) direction 
in Lai.85Sr.i5Cu04 would determine whether the gap is predominately s-wave or 
d-wave. We also propose an experiment to determine the gap at each point on a 
quasi-two-dimensional Fermi surface. 
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The anisotropy of the energy gap at the Fermi surface contains clues to the pairing 
mechanism operating in a superconductor. An s-wave gap^ successfully describes phonon- 
mediated pairing, which causes superconductivity in ordinary metals. A p-wave gap^ 
describes spin-flip-mediated pairing in helium-3. Current work suggests that the pairing 
in heavy fermion materials is mediated by anti-ferromagnetic fluctuations and is predom- 
inately d-wave^ . Theoretical work on high-Tg materials has produced theories too diverse 
and numerous to list, even as classes, but all must address the form of the gap function. 
Therefore we expect that accurate measurements of the gap would provide a strong test 
which any theory of high-Tg materials must pass. 

The experiments available to look for gap anisotropy provide limited information. 
These experiments can be divided into two classes: those which provide some average of 
the gap magnitude over all or a large part of the (normal) Fermi surface, and those which 
identify nodes in the gap on the Fermi surface. 

Reflectivity^, tunneling current^, ultrasonic attenuation^ and thermal conductivity^ 
measurements, as well as a variety of other transport measurements, yield some average 
value of the gap magnitude. In addition to their fundamental limitation - providing an 
average over the gap - there are other difficulties characteristic of each technique, which 
are discussed in the above references. Identifying nodes on the Fermi surface is done 
by measuring the temperature dependence of quantities such as ultrasonic attenuation^, 
specific heat^, magnetic penetration depth^^ or magnetic relaxation rate^^. Unfortunately, 
the exponential temperature dependence, that would rule out nodes, and the power-law 
dependence, that would indicate nodes, are difficult to distinguish. And even if the evidence 
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for a node is convincing, its focation is hard to determine because of the previously- 
mentioned averaging effects. 

We present a new method of analyzing measurements of phonon dispersion curves 
which should determine whether the gap in a high-Tc material is well described by a 
predominately s-wave or d-wave function. This method should also locate any nodes 
on the Fermi surface. In addition, recognizing that gaps which are approximately s, p, 
or d-wave are special cases of an anisotropic gap, we propose an experiment which can 
determine the gap magnitude everywhere on the Fermi surface. Both of these techniques 
require accurate measurements of phonon dispersion curves, which can now be made with 
inelastic neutron scattering^^. 

The gap at the Fermi surface and the Fermi surface shape completely determine the 
minimum energy required for a phonon to decay into two Bogoliubov quasiparticles. For 
phonon wavelengths smaller than the coherence length, the minimum-energy excitation 
is the creation of both quasiparticles on the Fermi surface. The positions on the Fermi 
surface where the quasiparticles are created can be found by placing the phonon momentum 
vector in the Brillouin zone so that its head and tail are on the Fermi surface (see Figure 
1). The quasiparticles are created at these points. The Fermi surface shape is described by 
a function kpiO), which is the vector whose tail is at the center of the Brillouin zone, whose 
head lies on the Fermi surface, and whose angle to a fixed axis is 0. Given the phonon 
momentum and the function kpiO), one can find the angles (p and (p where q's head and 
tail lie. Therefore we have the following equation, which should be viewed as determining 
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(j) and ^, given q and the function kpiO)' 

q^kF{4>)-kF{^). (1) 

The placement of the phonon momentum, and thus (j) and 0, are usually two-fold degenerate 
for a two-dimensional Fermi surface such as La2-xSrxCu04's. If the Fermi surface and the 
gap magnitude are inversion symmetric, this degeneracy is trivial and can be considered 
defined only from to tt. La2_xSrxCu04's Fermi surface is inversion-symmetric, and an 
inversion-symmetric gap magnitude conforms to current expectations-^^. We will discuss 
situations with non-trivial degeneracies at the end of this Letter. 

The minimum, or threshold, energy is the sum of the gap magnitude at the two points 
A;i?(0) and on the Fermi surface. These two points are uniquely determined by q and 

the function kpiO), so the threshold energy, denoted A, can be considered as a function 
of q and a functional oi kp- In the remainder of this Letter any dependence ori kp will 
be considered implicit. This paragraph, therefore, can be summarized by the following 
equation: 

A(g-) = |A(</.)| + |A(0)| (2) 

where A(0) is the gap at the point on the Fermi surface where the head of A;_p(0) is. As will 
be discussed later, A depends on q in markedly different ways if the gap is d-like instead 
of s-like. 

We will now present a way to determine the shape of the surface defined by A(g) in 
energy-momentum space. Phonons with energies below the surface, which we will call the 
threshold surface, should in general live longer than those above it. However, a variety of 
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other influences on the phonon hfetimes blur this distinction. Fortunately, near and just 
above the threshold surface, the effect on the phonon lifetime due to the newly opened 
quasiparticle channel becomes strongly amplified. The lifetimes and frequencies of phonons 
are usually plotted as a function of momentum magnitude, q, with the direction q fixed. 
With q fixed, A(g) defines a threshold line instead of a surface. Our proposal to determine 
is to observe anomalous behavior in the lifetimes and frequencies of the phonon 
branches which cross the threshold line. 

The transition from below to above the threshold line is sharp. Consider for illustrative 
purposes a weak-coupling BCS model where other influences on the phonon lifetime are 
approximated by a constant lifetime Tq. More realistic models will be discussed at the 
end of this Letter. We will follow a phonon dispersion curve which crosses the threshold 
line, beginning with a qi such that hijj{qi) < A (51) and concluding with a qs such that 
?iu>{qs) > ^{qs)- At the threshold momentum q2, the phonon frequency is on the threshold 
line: fiuj{q2) = K{q2). This model predicts the lifetime is Tq for qi < q < q2. The lifetime at 
q = q2 will vanish. At and slightly above the threshold line, the lifetime will be proportional 
to — (^2)^^^ and eventually will approach Tq if qs is sufficiently larger than q2. 

This anomalous behavior suggests that points on the threshold surface can be identi- 
fied by measuring the lifetimes of phonons and determining where in q, co space they drop 
precipitously. These points satisfy the equation 

M0 = |A(0)| + |A(0)| = A(g). (3) 

Equations (1) and (3) will be referred to as the kinematic equations. Each time a phonon 
dispersion curve crosses the threshold surface, one can determine a point on the surface. 
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Since there are a finite number of dispersion curves, A{q) will be measureable at a few 
points. This does not determine the function A(^) completely, but suffices to distinguish 
predominately s-wave, p-wave, or d-wave superconductivity. 

Figure 2 is a plot of the threshold line A{q) in the (100) direction for a modified 
tight-binding Fermi surface^^ (shown in Figure 1) with electron occupation 0.85. This is 
an approximation to the Fermi surface of Lai.85Sr.i5Cu04, which has the highest Tc for 
La2-xSrxCu04. Plotted is the threshold line for a d-wave gap 

A{e) = Ao [cos {kF{e)cose) - COS {kF{e)sme)] (4) 

and an s-wave gap A(^) = Aq. We draw the reader's attention to the node in A{q) near 
the half-way point to the zone boundary. The only way to have a node in the threshold 
line is to have that phonon momentum connect two nodes on the Fermi surface (shown as 
on Figure 1). So, this node in the threshold surface is a signature of nodes in the gap 
function. 

The number of (100) phonon branches which cross the threshold in the d-wave case is 
approximately equal to the number of zone center phonons with energies less than twice 
the maximum gap. High-Tc superconductors, therefore, are good materials to examine for 
this signature because of their large gap. 

The gap function magnitude can be determined everywhere on the Fermi surface 
by performing a different experiment. For a particular phonon branch, the kinematic 
equations have solutions which form curves in q,uj space. Assume that, given all these 
curves, at least one solution to the kinematic equations exists for each value of ^. It requires 
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quite a pathological case for this not to be true. For each pick one solution and identify 
it with (j). Knowing a solution to the kinematic equations means one knows the values of 
the quantities oj{q) and (f). So these quantities can now be considered functions of 0. Define 
the Fourier coefficients of a function of (f) as follows: /(0) = Ylm hnCOs{n(j)) + /2nSin(n(?!)). 
Equation (3), the energy equation, becomes the following set of coupled linear equations 
for the Fourier coefficients |A|j„ of the gap magnitude: 



and the other Cinjm involve the other three combinations of trigonometric functions. The 
inversion symmetry of the gap magnitude forces |A|^^ = for n odd. Matrix inversion of 
Equation (5) determines the |A|i^, and thus the gap magnitude everywhere on the Fermi 
surface. 

For a two-dimensional Fermi surface that is not inversion symmetric, or that is not 
convex, there are sometimes two or more non-trivial solutions to Equation (1) for q. The 
kinematic equations no longer always uniquely determine (p and in terms of q, therefore, 
the threshold energy becomes the minimum of |A(0)| + |A(0)| for all possible 0,0 
pairs. This newly defined threshold surface will still have structure if the gap is anisotropic. 
In particular, nodes in the gap function will be readily visible, with the threshold surface 
plunging to zero for a momentum connecting two nodes. The inversion procedure which 
culminates in Equation (5) remains unchanged if one is assured that any q, u solution to 
the kinematic equations which is used in the procedure only has a trivial degeneracy. 




(5) 



jm 



where 




(6) 
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For a three-dimensional Fermi surface this inversion procedure is no longer possible. 
In this case for a momentum g, Equation (2) has an infinite number of solutions for and 
(j), which are now solid angles. A direct algorithm for evaluating |A(6')| using kpiO) and 
the solutions to the kinematic equations is no longer possible because a fit of |A|j„ now 
requires multiparameter root finding. However the concept of a threshold surface remains 
valid and the nodes of it still indicate a momentum which joins gap nodes on the Fermi 
surface. 

Only a simple BCS model has been presented, but a calculation including strong 
electron-phonon coupling and Coulomb interactions^^'^^ indicates that the lifetime retains 
large anomalous features. In this calculation, the lifetime at g = ^2 does not reach zero 
and the discontinuity at this threshold momentum is smoothed. The rise of the phonon 
lifetime for q > q2, however, is much faster than in the weak-coupling case. Because of 
these competing efi^ects it is difficult to determine whether the transition from below to 
above the threshold surface would be easier or harder to observe in this more realistic 
model than in the weak-coupling model. 

The coherence factor in the electromagnetic response function provides a minor com- 
plication. For a system with the d-wave gap of Equation (4), this factor vanishes for a 
phonon whose wavevector is parallel to the (110) direction and whose head and tail lie on 
the Fermi surface. For such a phonon, therefore, the minimum-energy excitation does not 
contribute to its self-energy. This phonon's lifetime would still decrease upon crossing the 
threshold line, but would not be discontinuous at the threshold. We do not expect this 
behavior could be observed, considering current experimental capabilities. The solution to 
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this problem is simply to consider threshold lines in other directions, as done in Figure 2. 
Only particular phonon directions are affected, therefore this complication does not affect 
the inversion procedure described above for the gap magnitude. 

Since an extremely sensitive probe of surface phonons exists in thermal energy inelastic 
helium scattering-*^^, we remark that similar arguments to those presented in this paper 
may lead to gap measurement techniques involving surface phonons. 
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Figure 1. This is the Fermi surface for a one-band model with nearest-neighbor and 
next-nearest-neighbor hopping, taken from Ref. 14. The dispersion relation is 

e{k,9) = —2t{cos{kcos9) + cos(A;sin^)) — 4t'(cos(/i;cos^)cos(A;sin^)) 

with t — 430meV and t' — — 70meV. The electron filling factor is 0.85. qn connects two 
nodes on the Fermi surface if the gap is of this d-wave form. 

Figure 2. This is a plot of the threshold line in the (100) direction for the Fermi 
surface depicted in Figure 1. The solid line is the line for the d-wave gap A(^) = 
Ao [cos(/cf(6')cos6') — cos(A;F(6')sin6')] and the dashed line is for the isotropic A(^) = Aq 
case. 
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